Entanglement as a signature of quantum chaos 



Xiaoguang Wang, 1 ' 2 Shohini Ghose, 3 Barry C Sanders, 2 ' 3 and Bambi Hu 1,4 

1 Department of Physics and Center for Nonlinear Studies, 
Hong Kong Baptist University, Hong Kong, China. 
2 Australian Centre of Excellence for Quantum Computer Technology, Department of Physics, 
Macquarie University, Sydney, New South Wales 2109, Australia. 
3 Quantum Information Science Group, Department of Physics and Astronomy, 
University of Calgary, Alberta T2N 1N4, Canada. 
^Department of Physics, University of Houston, Houston, Texas 577204-5005, USA. 

(Dated: February 1, 2008) 

We explore the dynamics of entanglement in classically chaotic systems by considering a multiqubit 
system that behaves collectively as a spin system obeying the dynamics of the quantum kicked top. 
In the classical limit, the kicked top exhibits both regular and chaotic dynamics depending on the 
strength of the chaoticity parameter k in the Hamiltonian. We show that the entanglement of 
the multiqubit system, considered for both bipartite and pairwise entanglement, yields a signature 
of quantum chaos. Whereas bipartite entanglement is enhanced in the chaotic region, pairwise 
entanglement is suppressed. Furthermore, we define a time-averaged entangling power and show 
that this entangling power changes markedly as k moves the system from being predominantly 
regular to being predominantly chaotic, thus sharply identifying the edge of chaos. When this 
entangling power is averaged over initial states, it yields a signature of global chaos. The qualitative 
behavior of this global entangling power is similar to that of the classical Lyapunov exponent. 

PACS numbers: 05.45.Mt, 03.65.Ud, 03.67.-a 



I. INTRODUCTION 

Quantization of classical chaotic systems has long been 
of interest because of the challenges in identifying quan- 
tum signatures of systems that, in a classical limit, ex- 
hibit chaotic behavior. Various signatures have been 
identified, such as the spectral properties of the generat- 
ing Hamiltonian p], phase space scarring 0, hypersen- 
sitivity to perturbation Q, and fidelity decay Q, which 
indicate an underlying chaotic presence in the quantum 
dynamics. Here, we investigate the issue of quantum sig- 
natures from the perspective of entanglement: as entan- 
glement is at the heart of quantum mechanics and a cru- 
cial resource for quantum information processing 0, Q , 
the entanglement inherent in quantum chaotic systems 
could provide a valuable approach to studying decoher- 
ence and quantum chaos U H, 0, El EH • Further- 
more, quantum chaos could be seen as an engine for gen- 
erating entanglement; indeed quantum chaos could lurk 
in quantum information processing [Til. IbH HEj as a dele- 
terious or perhaps even as an advantageous effect. We 
study entanglement by considering a finite multipartite 
system, whose collective dynamics obey chaotic Hamil- 
tonian dynamics in the classical limit. 

Previous studies of entanglement in chaotic systems 
0EH0[H[ll[ll0 have explored bipartite entan- 
glement in pure states, entanglement of qubits in a mul- 
tiqubit system and average entanglement or entangling 
power. In this study we explore these different types of 
entanglement in a single simple system, namely the quan- 
tum kicked top (QKT) IH HE Ell. This enables 
us to compare the degree to which these different types 
of entanglement display signatures of chaos. 



An advantage of the dynamics of the QKT is that it 
obeys a spin algebra symmetry. This spin system can 
thus be regarded as a composite of distinct spin-half par- 
ticles thereby admitting a multiqubit interpretation. Our 
system thus allows us to study and compare pairwise en- 
tanglement between two qubits as well as bipartite en- 
tanglement between the two qubits and the rest of the 
qubits. The Hilbert space for the QKT is finite and the 
Poincare section of the phase space compact, allowing 
analyses of quantum and classical dynamics uncompli- 
cated by truncation issues. The QKT is well studied and 
understood thereby simplifying the analysis of the role 
of entanglement in the system. Finally, the QKT pos- 
sesses a parity symmetry that allows coherent quantum 
tunnelling to occur for states localized at classical fixed 
points 20] ■ 

The Hamiltonian evolution may increase the entangle- 
ment of the multipartite system, initiated in a collective 
spin coherent state 21]. For this analysis it suffices to 
employ two entanglement measures. For bipartite entan- 
glement, where the multipartite system is divided into 
two subsystems, entropy of a subsystem is used to quan- 
tify the degree of entanglement between the two subsys- 
tems. Pairwise entanglement, on the other hand, consid- 
ers the degree of entanglement between two qubits traced 
over all remaining qubits and is quantified by the concur- 
rence [2l|23|. 

We present general results for both bipartite and pair- 
wise entanglement in the multipartite QKT and demon- 
strate that these entanglement measures reveal strong 
signatures of the classical chaos features corresponding 
to the onset of chaos and to the edge of chaos [24J , which 
is the boundary between regular and chaotic regimes in 
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the classical phase space. We have studied the behav- 
ior of the linear entropy and the concurrence for specific 
initial states as well as the dynamics of these quantities 
when averaged over all initial states. We show that the 
linear entropy increases more rapidly for an initial state 
centered in a chaotic region of the classical phase space 
than one centered on a fixed point. This agrees with the 
behavior of the linear entropy observed in other chaotic 
systems, supporting the conjecture that classical chaos 
can enhance bipartite entanglement 0, @ ■ Furthermore, 
we show that the pairwise entanglement as measured by 
the concurrence also reveals a dramatic change for a spin 
coherent state whose mean traverses the edge of chaos 
on its transit through chaotic and regular regions of the 
phase space. Contrary to the linear entropy, the concur- 
rence rapidly decreases for an initial state located in the 
chaotic region. 

While the linear entropy and concurrence can reveal 
the local chaotic and regular structures in phase space, 
the entangling power, which is the averaged bipartite or 
pairwise entanglement, can be used to identify the edge 
of chaos and quantify the onset of global chaos, much 
like the classical Lyapunov exponent. We show that the 
entangling power greatly increases as the chaoticity pa- 
rameter k is increased, and the corresponding classical 
kicked top makes the transition from predominantly reg- 
ular to predominantly chaotic behavior. In particular, 
the behavior of the average linear entropy is qualitatively 
similar to that of the classical Lyapunov exponent, thus 
revealing a signature of a global feature of the classical 
chaos. 

The paper is organized as follows. In Sec. II, we in- 
troduce the QKT and its classical dynamics, and intro- 
duce bipartite and pairwise entanglement measures. In 
Sec. Ill, we study in detail dynamical evolutions of bipar- 
tite and pairwise entanglement, and examine the edge of 
quantum chaos, the onset of quantum chaos via the en- 
tangling power. We conclude in Sec. IV. 

II. BACKGROUND 
A. Quantum kicked top 

The QKT is described by the Hamiltonian [Tlll7lll^| 



H=—J 2 z +pJ y £ S(t-nr), (1) 

n— — oo 

where J a (a £ {x,y,z}) are spin operators and states 
are restricted to irrep j for which J 2 = j(j + 1). r is 
the duration between periodic kicks, p is the strength of 
each kick (which is manifested as a turn by angle p), and 
k is the strength of the twist. The Hamiltonian is an 
alternative sequence of twists (J 2 term) and turns (J y 
term). The QKT describes a spin system, which can be 
comprised of multiple systems of lower spins. For {(7i a } 



the Pauli operators for the i qubit, a collective spin 
operator, satisfying the usual su(2) algebra, is given by 

JV 
i=l 

An example of using multiple qubits to simulate the QKT 
has been presented for trapped ions [2S|. 

A standard dynamical description of the QKT is via 
the Floquet operator 

F = exp ^^—Jl^j ex P {-VpJy) , (3) 

where the energy is rescaled so that r = 1, and p = n/2 
are henceforth assumed. The orthogonal eigenstates of 
F, denoted by {|$ m ) : — j < m < j}, which satisfy 

F\$ m ) =exp(i$ m )|$ ro ), (4) 

provide a convenient basis for stroboscopic evolution. An 
arbitrary state |^(0)) evolves to 

j 

|*(n)> = F n |*(0)) = *m(0)exp(in$ m )|$ ro ). (5) 

with * m (0) =1*jJ*(0)>. 

The QKT [Tq, [13 is chaotic in the classical limit. For 
integrable systems, it is well known from semi-classical 
theory that the classical actions can be associated with 
corresponding regular eigenstates of the quantum system 
with a well defined quantum number. This correspon- 
dence breaks down in chaotic systems |26| . In quasi- 
integrable systems with a mixed classical phase space 
of regular and chaotic regions, some of the eigenstates 
can still be associated with local actions in the regular 
regions with corresponding discrete eigenenergies. The 
remaining eigenstates result in an irregular energy spec- 
trum corresponding to the chaotic region [23. We show 
here that this underlying regular and chaotic energy spec- 
trum of the Floquet eigenstates of the QKT is reflected in 
the dynamics of the entanglement, depending on whether 
the initial state is in the regular or chaotic region of the 
classical phase space. 

The classical limit of the QKT is obtained by express- 
ing X = (J x /j) and similar for Y and Z and factorizing 
all moments such as (J x Jy/j 2 ) = XY to products of first- 
order moments. Then the classical equations of motion, 
obtained from the Heisenberg operator equations of mo- 
tion and applying the factorization rule above, are given 

by na 

X' = Zcos(kA) + y"sin(«;A), 

Y' = -Zsin(«;A) +Fcos(kA), (6) 

Z' = -X 

The stroboscopic evolution described by Eq. @) can be 
represented in a phase space given by a sphere S 2 of unit 



3 




FIG. 1: The stroboscopic phase space dynamics of the classi- 
cal kicked top for re = 3. Three hundred stroboscopic trajec- 
tories are plotted, each for a duration of 300 kicks. 



radius. The classical, normalized angular momentum 
variables (X, Y, Z) can be parametrized in polar coordi- 
nates as (X,Y,Z) = (sin cos sin sin 0, cos 8), where 
and 4> are the polar and azimuthal angles, respectively. 
Thus the mapping domain is essentially two-dimensional. 

The stroboscopic dynamics of the classical map is 
shown in Fig. 1. In the plot, we choose the chaoticity 
parameter n — 3 which yields a mixture of regular and 
chaotic areas of significant size. Elliptic fixed points sur- 
rounded by the chaotic sea are evident. Two such elliptic 
fixed points have coordinates (9, (f) = (2.25, —2.51) and 
(6, <f) = (2.25,0.63). As we will see, this phase space 
structure of the classical kicked top determines behavior 
of quantum entanglement in the QKT. 



B. Entanglement measures 

Pure-state bipartite entanglement has been calcu- 
lated in previous studies that explore connections be- 
tween quantum entanglement and underlying classical 
chaos [H, Ig L flCX ll ll [l2( • Quite recently, Bettelli and She- 
pelyansky jjjj studied the behavior of the concurrence in 
a system exhibiting quantum chaos, and found that the 
underlying classical chaos leads to exponential decrease 
of the concurrence down to some residual values. This 
result shows that the concurrence is very sensitive to the 
onset of chaos, and understanding entanglement of these 
systems may help to control quantum chaos and suppress 
its negative effects in quantum information processing. 

We consider bipartite entanglement between a pair of 
qubits selected from a symmetric multi-qubit state and 
the rest of the system, as well as pairwise entanglement 
between the two qubits. Once we obtain the two-qubit 
reduced density matrix, the entanglement can be readily 
calculated. By expressing the reduced density matrix p\2 
for qubit 1 and 2 in terms of the expectation values of the 
collective operators, all elements of p\2 are conveniently 



obtained HJ. 

Our system governed by the QKT Hamiltonian is a 
composite system, and remains in a pure state at all times 
if we initially choose a pure state. For pure states, bipar- 
tite entanglement is well-defined and can be quantified 
by entropies of either subsystem. For convenience, we 
adopt the linear entropy as the entanglement measure, 
which is defined as 

E = l-Trx(pl), (7) 

where p\ is the reduced density matrix for the first sub- 
system. The maximum linear entropy for a pure state of 
a bipartite d\ x d,2 system, is given by 1 — 1/ min(<ii, c^)- 
While we may choose other entropies such as the von 
Neumann entropy as our entanglement measure, the 
qualitative results are in general independent of choice 
of entropies for pure states. Moreover, the linear entropy 
and the von Neumann entropy are two limiting cases of 
the Renyi entropy |2j| , they are thus interrelated and one 
can be used to estimate the other [3(], |^ . 

Given our iV-qubit system, we consider another type 
of entanglement, the pairwise entanglement, i.e., the en- 
tanglement between a pair of qubits. When TV > 3, the 
pair of qubits can be in a mixed state. Entanglement 
for a mixed state p\2 is quantified by the entanglement 
of formation. Specifically, for a pair of qubits, entangle- 
ment is equivalent to the n on-p ositivity of the partially 
transposed density matrix [32| . Alternatively, one can 
use the concurrence |22l l23j to quantify the pairwise en- 
tanglement. The concurrence is defined as 

C = max {0, Ai - A 2 - A 3 - A 4 } , (8) 

with the quantities Ai being the square roots of the 
eigenvalues in descending order of the matrix product 
Pi2(ciy 8> <J 2y)Pi2(< 7 iy ® &2 y )- P12 denotes the complex 
conjugate of p\2- The value of the concurrence ranges 
from zero for an unentangled state to unity for a maxi- 
mally entangled state. 

III. ENTANGLEMENT AND QUANTUM 
CHAOS 

We present here our studies of the entanglement dy- 
namics of our 7V-qubit system governed by the QKT Q , 
with the relevant angular momentum operators J a given 
by the collective operators. If we choose the initial pure 
state to be symmetric under exchange of any qubits, 
then the state vector at any later time is also symmet- 
ric. Thus, we can describe the state of the A^-qubit 
system in terms of the orthonormal basis { \j, to) ; (to = 
-j,-j + 1,. with j = N/2. The states {\j,m)} 
are the usual symmetric Dicke states |33j. State \j,—j) 
is not entan gled , whereas state \j, —j + 1), the so-called 
W state |3~H l35j| . is pairwise entangled with concurrence 
C = 2/N'. 

To connect the quantum and classical dynamics of the 
kicked top, we choose the initial state to be the spin 
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FIG. 2: Dynamical evolution of the linear entropy for initial 
SCS \0 = 2.25, d> for ft = 3 and N = 50. 



FIG. 3: Dynamical evolution of the concurrence. The pa- 
rameters are the same as those of Fig. 2. 



coherent state (SCS) \\9,(j>) = R(9,4>)\j,j);-ir < <j> < 
tt,0 < 9 < tt}, with 

R(8,(j)) = exp{i9[J x sin <j> — J y cos </>]}. (9) 

Mean of J / j is 

(0, </>| J/j \9, 4>) = (s'm9 cos </>, sinf? sin0, cos 0). (10) 

The initial SCS can be rewritten as a multi-qubit prod- 
uct state, and thus exhibits no entanglement (zero linear 
entropy and concurrence). 

A. Dynamics of entanglement 

We start by exploring the dynamics of entanglement 
for initial states with a mean value in four different re- 
gions of the phase space, specifically a fixed point, an in- 
tegrable (or KAM) region, a chaotic region, and the bor- 
der between the integrable and the chaotic region known 
as 'the edge of chaos' [24] . We are also interested in the 
behavior as the chaoticity parameter n is varied. We start 
with the choice k = 3, which exhibits large integrable and 
large chaotic regions, and we select four states localized 
in the four regions mentioned earlier. 

For convenience we fix 9 = 2.25 and vary </>: this 'line of 
latitude' on S 2 includes all four regions we are exploring. 
An elliptic fixed point arises at <j) = 0.63, a point in the 
regular region occurs at (f> = 0.90, one edge of chaos can 
be seen at (f> = 1.05, and a point well in the chaotic sea 
is located at (f> = 2.00. The states with means at each 
of these (9, <fi) points in the phase space are chosen to be 
SCSs. These states are minimum uncertainty states and 
are well localized around the four chosen points in phase 
space. 

We numerically compute bipartite entanglement be- 
tween two qubits and the other N — 2 qubits, and dis- 
play the numierical results in Fig. 2 as linear entropy 
versus time increasing parameter n. We observe that the 
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FIG. 4: Long time behavior of the concurrence for the regular 
case (9 — 2.25, (f> = 0.63). The other parameters are the same 
as those of Fig. 2. 



entanglement is enhanced for the initial state centered 
in the chaotic region after a short time. Initially, the 
linear entropy is zero, and as the dynamics evolve, the 
entropy increases slowly when the wave packet is cen- 
tered in the regular region, whereas it exhibits a rapid 
rise when centered in the chaotic region. The curve with 
<f> = 1.05 displays the intermediate behavior. Further- 
more, the entanglement for the state initially with mean 
at a fixed point displays a periodic modulation that is 
absent in the evolution of the entanglement for the state 
initiated in the chaotic region. This periodic modula- 
tion is an indicator of the underlying regular classical 
dynamics and corresponding regular energy spectrum of 
the Floquet eigenstates [22| ■ 

Figure 3 shows the dynamical behavior of the concur- 
rence. Just as in the linear entropy case, we see a rapid 
change in the concurrence for a state initially centered 
in the chaotic sea. However, our numerical results sug- 
gest that the initial state centered in the fixed point will 



lead to large pairwise entanglement production, which is 
opposite to the case of the bipartite pure-state entangle- 
ment production, in which the classical chaos enhances 
the production of bipartite entanglement. Of particu- 
lar interest are the collapses and revivals in the evolu- 
tion of the concurrence for a state centered on the fixed 
point. Figure 3 (left-up one) shows a revival at n = 52. 
Additional revivals occur at n = 113, n = 183 and so 
on. For the chaotic case, we cannot observe the revival 
phenomenon. This quasiperiodic behavior in the regu- 
lar region indicates that the SCS in the regular region 
has a finite support over the basis set of regular Floquct 
eigenstates. 

In Fig. 4 we plot the long-time behaviors of the con- 
currence for the regular case, and observe multiple col- 
lapse and revivals. At long times the revivals become 
sparse, and finally the concurrence reduces to zero. We 
see that the concurrence is very sensitive to quantum 
chaos, which is consistent with the observations of Bet- 
telli and Shepelyansky [13j in their studies of concurrence 
between qubits during the operation of an efficient multi- 
qubit quantum algorithm. However, unlike their system 
in which the concurrence reached a finite residual value, 
in our QKT model, the concurrence disappears at very 
long times; the difference arises because we assumes sym- 
metrised multi-qubit states, whereas Bettelli and Shep- 
elyansky allow this symmetry to be broken. 



B. Edge of quantum chaos 

The edge of chaos [24J is an important issue in the 
study of quantum chaos. In classical chaos, the edge of 
chaos is a fractal boundary separating the regular and 
chaotic regions. However, this fine-grained fractal struc- 
ture does not translate well into the quantum domain. 
Recently, it was found that the edge of chaos is charac- 
terized by a power law decrease in the overlap between a 
state evolved under chaotic dynamics and the same state 
evolved under a slightly perturbed dynamics (24|. Here, 
we study the edge of quantum chaos from the perspective 
of entangling powers, which are defined to be either the 
maximal or the mean entanglement that the evolution 
operator can generate over all initial states [3(| |3]| • Al- 
ternatively, given a fixed initial state, we may ask what is 
the maximal and the mean entanglement that the opera- 
tor can generate over all time. In general state-averaging 
and time-averaging are inequivalent and so the two meth- 
ods yield different results. 

In strongly chaotic systems, the two definitions con- 
verge due to nearly ergodic dynamics. In this study, we 
explore both the entanglement averaged over all time as 
well as the entanglement averaged over initial states. We 
begin our analysis with the average over time. In prac- 
tice, for numerical purposes we consider a finite time do- 
main. We define time-averaged entanglement power as 
the average linear entropy or average concurrence over a 
time interval T (which should be much longer than other 
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FIG. 5: Mean linear entropy and mean concurrence against 
azimuthal angle 4>. The mean linear entropy is plotted for 
N = 20 (solid line), N = 50 (dashed line), and N = 100 
(dotted line); the mean concurrence is plotted for N = 100 
(solid line), N = 200 (dashed line), and N = 300 (dotted 
line). Parameters 6 and k are set to 8 = 2.25 and k = 3. 
The time average is over 200 steps, and the subplots below 
are partly enlarged versions of above ones. 



time scales) as follows: 



E T = i jTdi E(t), C T = 1 jTdt C(t). (11 



For local unitary operations, the above quantities are nec- 
essarily zero. 

We fix the polar angle 8 = 2.25 of the SCS as before, 
and vary the azimuthal angle <f> from — tt to tt. The center 
of the SCS wave packet thus commences in the chaotic 
region and passes through two regular islands. Figure 5 
displays the time-averaged mean linear entropy i/rand 
mean concurrence Ct as a function of the azimuthal an- 
gle <f>. When the azimuthal angle goes from — tt to the 
first regular region, the linear entropy decreases until it 
reaches a minimum which approximately corresponds to 
the fixed point (6,<f>) — (2.25,-2.51). Subsequently the 
mean entropy increases to a flat larger area correspond- 
ing to the chaotic region. 

In contrast to the behavior of the mean linear entropy, 
the mean concurrence reaches a maximum approximately 
at the fixed point. When (f> increases from — tt, the mean 
concurrence first decreases slowly, and then exhibits an 
abrupt increase to a maximal value. Thos turning point 
sharply displays the edge of chaos. Another turning point 
is obvious from the figure. 

We also observe that the larger the number of qubits, 
the wider the regular region. In the large N limit, the 
regular region will coincide with that of corresponding 
classical chaos of Fig. 1. 

We further calculate the mean linear entropy and the 
mean concurrence as a function of 9 and <p. The contour 
plots are shown in Figs. 4 and 5. Comparing Figs. 1 and 
4, we observe that these two figures closely match each 



FIG. 6: Contour plot of mean linear entropy against cf> and 
6. The mean is over 200 kicks and other parameters are the 
same as those in Fig. 2. 
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FIG. 7: Contour plot of mean concurrence against <f> and 9. 
The average is over 200 kicks and other parameters are the 
same as those in Fig. 2. 



other. Specifically, the four islands of Fig. 4 are evident, 
reflecting the four stable islands in the classical phase 
space. Comparing Figs. 1 and 5, the four stable islands 
of Fig. 1 closely match those of Fig. 5. Thus, we have a 
good classical-quantum correspondence. 



C. Onset of Chaos 

In the previous section, the time-averaged entangle- 
ment revealed clearly whether the initial state was in the 
regular or chaotic region. Here, we are concerned with 
global properties of a chaotic system. For classical sys- 
tems, the onset of global chaos can be quantified by cal- 
culating the global Lyapunov exponent. Here we define 
the following entangling power to quantify the onset of 



FIG. 8: Entangling power er (square line), ct (diamond line), 
and the Lyapunov exponent A (circle line)agaist k. The plot 
of Lyapunov exponent corresponds to fig. 1 of Ref. |3q| . The 
parameter N = 36 and time average is over 50 kicks. 

global quantum chaos: 

e T = i J^J dtE(t,0,if>), 

c T = ±Jdv^ dtC(t,0,<f>), (12) 

where d/i = dfi(9, </>) is the Haar measure. Like the global 
classical Lyapunov exponent, this state-averaged entan- 
gling power characterizes global properties of the QKT. 
The entangling power characterizes the entangling capa- 
bility of the K-dependent Floquet operator. 

Figure 8 shows the entangling powers and the Lya- 
punov exponent versus parameter n. We observe that 
when k — kq ~ 2.4, the entangling power ex exhibits a 
rapid increase and saturates beyond k — k\ ps 5. The 
rapid increase of the entangling power signifies the on- 
set of quantum chaos, and the saturation implies that 
global chaos has occurred. Between Kg and K\, when 
there are still regions of regular islands in the chaotic 
phase space, a mixture of regular and chaotic behavior is 
expected. In contract to ex, ct becomes very small for 
n > kq, which is also an indicator of onset of quantum 
chaos. Note that Ct has a peak which results from the 
competition between the entangling power of the QKT 
Floquet operator and the inherent chaos. On one hand, 
increasing k will enhance the entangling power, and on 
the other hand, the inherent quantum chaos suppresses 
the pairwise entangling power, thus leading to the peak. 
For ex, increase of k and the quantum chaos both en- 
hance the linear entropy, and thus no competition exists 
and no peak appears. 

IV. CONCLUSIONS 

We have investigated a multi-qubit system whose col- 
lective Hamiltonian dynamics are chaotic in the classical 
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limit. We studied the particular example of the quantum 
kicked top, which is a well-studied example of quantum 
chaos with the advantages of having a finite-dimensional 
Hilbert space (thereby obviating the need for truncation 
that arises in infinite-dimensional Hilbert spaces) and of 
involving only spin operators of no more than quadratic 
order. 

Although the collective dynamics are well understood, 
the underlying entanglement of the qubits that collec- 
tively make up the quantum kicked top is only just begin- 
ning to be understood. Here we have developed methods 
for studying the quantum kicked top, and these methods 
are applicable to more general systems. We have iden- 
tified bipartite and pairwise entanglement as two quite 
distinct measures to determine the entanglement in the 
system, and we have related the dynamics of these mea- 
sures of entanglement to chaotic features of the quantum 
kicked top in the classical limit; as examples, we have 
connected the features of the entanglement evolution to 
local properties such as whether the state is supported 
predominantly in the regular or chaotic region and also to 
global properties such as showing that entangling power 
averaged over states grows similarly to the global Lya- 
punov exponent growth for the classical chaotic system. 

We have assumed symmetric multi-qubit states 
throughout, and the entanglement properties studied 
here reflect this assumption. If the symmetrization con- 
dition is broken, different dynamics can be expected. For 
example, Bettelli and Shepelyanski ^| show a concur- 
rence that reaches a residual steady-state value. They 
explain this non-zero residue as being a result of sym- 
metry breaking in their system. In contrast, our system 



exhibits a decay of concurrence to zero. The assump- 
tion of symmetric states implies indistinguishability of 
the qubits. Thus even though entanglement may exist 
in the system, it may not be accessible as a useful tool 
for quantum information processing due to inherent in- 
ability to distinguish between the qubits. An alternative 
measure of entanglement could be an operational mea- 
sure that takes into account physical restrictions on ac- 
cessibility of the entanglement due to symmetries of the 
system |3!|. 

In summary, our work highlights the connection be- 
tween entanglement of a multi-qubit state whose collec- 
tive dynamics is chaotic in the classical limit and in- 
troduces valuable methods and measures for studying 
this entanglement. It would be worthwhile to investigate 
other quantum chaotic systems using concepts of time- 
averaging and global entangling power. Furthermore, it 
would interesting to define and compare systems in which 
the symmetrization results hold to systems where this 
symmetrization condition is broken. 
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